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Abstract. An important aspect of astrophysical MHD turbulence research is develop- 
ing diagnostics to connect simulations with the observable universe. Turbulent systems 
are by definition structurally complex in all fluid variables (density, velocity, and mag- 
netic field), such that they must be described statisticaUy. By developing and applying 
diagnostic tools to simulation data, it is possible to interpret empirical laws for the sta- 
tistical properties of observed systems in terms of fundamental dynamical processes, 
and to identify and calibrate robust probes of physical parameters that cannot be mea- 
sured directly. Using several different examples, I describe how structural diagnostic 
analyses have already yielded significant insights into the nature of turbulent molecular 
clouds. I review results from several different groups, and discuss directions for future 
diagnostics to enhance our understanding of cloud structure and constrain models of 
the evolutionary course that governs star formation. 



1 Introduction 

As the number, range, and depth of the papers in this volume witnesses, re- 
cent progress in modeling and understanding astrophysical MHD turbulence is 
impressive. Even with the intensive research of several groups over the past 
few years, however, many aspects of the fundamental turbulence phenomenon 
are not yet wholly understood - which makes for continuing excitement in this 
emerging discipline. In using the results of MHD simulations to interpret the 
dynamics of the interstellar medium, the technical challenges involved in numer- 
ically modeling and characterizing turbulence are compounded by astrophysical 
uncertainties in posing the numerical problem to be solved. For molecular clouds, 
open astrophysical questions include: 

• What is the source (original, and potentially, maintaining) of turbulence? 

• What is the mean magnetic field strength, and variation of mass-to-flux ratio, 
in molecular clouds? 

• What is the range of sizes and masses of molecular clouds in the Milky Way 
and other spiral galaxies? 

• How are clouds formed? how long do they live? how are they destroyed? 

From the point of view of defining an idealized problem for an MHD simulation, 
these astrophysical questions translate to uncertainties in the input spectral 
form (in space and time) of the turbulent driving, the value and variation of the 
plasma P parameter, the importance of self-gravity, and the initial and boundary 
conditions for the simulation. 
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The complexity of the turbulence phenomenon demands detail and variety in 
the analytical methods used to characterize its structure. For application to un- 
derstanding astronomical systems - where the physical inputs are uncertain, and 
only projected distributions are available - the eventual aim is to develop a set 
of simple, robust diagnostics of MHD turbulence that have direct connections to 
observable quantities. Potentially, there are many different avenues for this sort 
of analysis, and extensive exploration is required to determine what directions 
are most productive. Because large-scale numerical simulations of turbulence un- 
der interstellar conditions are only now becoming computationally practical, the 
present tasks include first the "forward" process of characterizing MHD turbu- 
lence obtained from new simulations with a range of parameter values, and then 
using these results to select and calibrate diagnostics for the "reverse" process 
of discriminating systemic parameters from observables. 

As examples of the process of developing diagnostics of molecular clouds' 
internal structure, kinematics, and magnetization from turbulent MHD simula- 
tions, my discussion here will focus on recent work analyzing and interpreting 
density and column density statistics (§2), properties and definitions of clumps 
(§3), linewidth-size relations (§4), and statistics of polarization maps (§5). I 
will update previous work (see also Q, 0), present several new 

results, and make connections to the conclusions of other groups. Chapters in 
the same volume covering topics most directly related to those addressed here 
include those by Crutcher, Heiles, & Troland; Nordlund & Padoan; MacLow; 
Cho, Lazarian & Vishniac; and Zweibel, Heitsch, & Fan. 



2 Density and Column Density Statistics 

In general, a turbulent velocity field leads to production of significant local den- 
sity variations in a compressible medium (i.e. a medium with Mach number 
M = I'rms/cs 3> 1, where Wi-ms is the turbulent velocity dispersion and Cg the 
sound speed). This is true regardless of the magnetic field strength, because in 
the case of a weak magnetic field (wa = Bq / ^/Anp Cg Vmis, where Bg is 
the mean magnetic field and p is the mean density), magnetic pressure forces 
are weak compared to ram pressure, and in the case of strong magnetic fields 
(va "frms), compression is unhindered along the mean field direction (and 
indeed enhanced by forces associated with gradients in B^, where B_l is the 
component perpendicular to Bq). 

If compression and rarefaction events are spatially and temporally inde- 
pendent, then for the case of an isothermal equation of state (approximately 
true under molecular cloud conditions), the resultant one-point density distri- 
bution function (often referred to as a "PDF" - probability distribution func- 
tion) is expected to obey a lognormal form (js^, |29|). If 1 -t- (5,; is the enhance- 
ment/decrement factor for density in the i*^ compression/rarefaction event af- 
fecting a given fluid element, then the density after N events will be 

P^pUl,il + 6.) (1) 
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SO that the logarithm of the density, 

log(p/p)=E^ilog(l + <5,), (2) 

will be the sum of independent random variables; by the Central Limit Theorem, 
this implies that \og{p/p) should obey a Gaussian distribution. 

The results of numerical simulations (3|],[0,|^) indeed bear out 

the expectation that a log-normal form for the volume density PDF prevails 
(at least away from the tails) under isothermal conditions. This result holds 
both for forced and decaying turbulence, and for simulations with varying mean 
magnetic fields. Figure 1 shows an example of the distributions of fractional 
volume and mass as a function of volume density for four forced-turbulence MHD 



simulations with = 5 (see [44 for details on the models), with comparisons 
to the lognormal functions with the same mean and dispersion. 

As the example in Fig. 1 shows, the average mass compression factor is 
relatively insensitive to the mean magnetic field strength. Although the minimum 
value of the mass-weighted mean (log(p/p))M increases (logarithmically) with 
the fast-magnetosonic Mach number My = Wrms/(Cs -I- w^)^/^, the scatter from 
"cosmic variance" is large enough that there does not appear to be a unique 
relation between A^f (or ; cf. |32| , |^ ) and {\og{p/p)) [Q. Typical values of 
the mean compression factor for matter in the simulations, {p/p)m, range over 3- 
6 for A4 = 5 — 9, consistent with the compression factor needed to excite the CO 
molecule - rendering molecular clouds observable - when the volume-weighted 
density is only n ~ 100 cm^'^ (e.g. |42| ). 

More directly observable than the distribution of volume densities in a cloud 
is its distribution of column densities, corresponding to line-of-sight integrations 
of density, S = J pds. Heuristically, one might expect the column density at any 
projected position to be determined by a series of independent (in space and 
time) compressions and rarefactions, similarly to the process of events described 
in (1). The difference for column density is that each event affects only a fraction 
/i < 1 of the line of sight, so that the column density is given by 

log(i;/^) =E,log(l + /.(5,) (3) 

instead of (2). The factor fi may be thought of as the ratio of the correlation 
length of a given compression/rarefaction event to the overall linear scale of the 
cloud along the line of sight. If the individual enhancement / decrement factors are 
independent random variables, then the resultant column density PDF should 
take on a log-normal form. Because each 1 -I- fiSi is closer to unity than 1 + Si, 
however, the mean and dispersion of \og(S / S) are expected to be smaller than 
the corresponding quantities for log{p/p). These expectations are indeed borne 
out by analyses of column densities in simulations, as shown by [^Tf ; distributions 
show a lognormal form (see also Fig. 2), and typical values of the mass- weighted 
mean column density are {S / S)m = 1-1 — 1-4 for = 5 — 9. 

From (3), note that if the factors fi are small (corresponding to having 
the dominant correlation length small compared to the size of the numerical 
box/physical cloud), then 5S/S = E/E — 1 Hifidi, implying a Gaussian 
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Fig. 1. Distributions of fractional volume (solid curves) and fractional mass (dashed 
curves) as a function of density for Macli-5 forced turbulence simulations with varying 
mean magnetic field strength characterized by /3 = c^/v\, for va = -Bo /\/47rp. Dotted 
curves show lognormal distributions for comparison. Each panel is also labeled with 
the mass- and volume-weighted mean compression magnitudes. 



distribution for 5S/S if the terms fi5i are independent random variables. Nu- 
merical evidence on how column PDFs transition from Gaussian to lognormal 
form as fi increases is presented in 

Preliminary comparisons between observed PDFs of column density in fil- 
amentary molecular clouds - obtained from stellar extinction data |^ - and 
simulated PDFs from turbulence models are very encouraging |^ . The data are 
consistent with lognormal distributions, with comparable width to those from 
simulations having turbulent Mach numbers and power spectra comparable to 
those in observed clouds. It remains to be seen how much more specific infor- 
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Fig. 2. Distributions of fractional area (histograms) as a function of column density 
for Mach-7 decaying turbulence, for two different projection directions, from l3 = 0.01 
simulation (see [Bl| for details). Solid curves show the corresponding lognormal distri- 
bution fits; dashed lines indicate exponential fits on high-column side, with slopes -1.1 
and -1 for the x and z projections. 



mation about a cloud can be learned from its column density. Although there 
were earlier some hopes that column density PDFs could help distinguish the 
mean magnetic field strength in clouds , models with matched Mach numbers 
"observed" from varying directions do not show strong or consistent correlations 
of the mean column density contrast with the value of /3 = c'^{Bq /4:np)~-^ ||3l|] . 
In spite of this insensitivity to the mean magnetic field strength, one might still 
hope to constrain the distribution of volume densities from the distribution of a 
cloud's column densities. The mathematical degeneracy between Si (the volume 
compression/rarefaction increment) and /,; (the spatial coherence length of an 
event) evident in (3) indicates that there is no simple inversion method. Methods 
that combine the column PDF with the spatial (two-point) correlation function 
in column density maps (reflecting fi, or more generally, the shape of the density 
power spectrum) may however be able to lift this degeneracy; this represents an 
important direction for future study. 

Because the population of the tails of the density and column density PDFs 
may be more affected by intermittency than the population near the peak, and 
because the equation of state may depart from isothermality in more overdense. 



optically thick regions (cf. |41 ,[g9[) observational departures from log- normality 
are more likely to occur there. Simulations (even with a uniform, isothermal 
equation of state) show a variety of behaviour in the tail distributions, although 
PDFs that are lognormal over more than three orders of magnitude are common 
from isothermal decaying-turbulence simulations. Fig. 2 shows, for example, that 
for column densities less than S/S = 4 (corresponding to some 98% of the 
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projected area of the particular simulation), a lognormal function is an excellent 
approximation. 

Over a limited range of column densities above the mean, lognormal PDFs 
can typically be fit by an exponential function as well (i.e. log[(iA/(iZ'] = Ci + 
C2S / E for Ci and C2 constants). Fig. 2 shows examples of local fits of this kind 
with slopes C2 = —1.1,-1. Similarities between this local exponential form in 
simulations and in molecular- line observations with emphasis on potential 
dependences of the slope on the largest density correlation scale, have been in- 
vestigated by 1^ (see also Q for related discussion of resolution effects) . It is 
not yet clear whether there is inherent physical significance in the local exponen- 
tial form, or whether it is primarily a convenient mathematical approximation 
to a lognormal on the high-column-density side of the distribution where most 
of the observable matter is found. 



3 Clumps in Turbulent Clouds 

A longstanding unsolved problem in astrophysics is what determines the stellar 
initial mass function (IMF). Many different physical processes could potentially 
affect the IMF; an abbreviated list includes: (i) turbulent stresses in a large- 
scale cloud producing non-self-gravitating clumps with a range of masses/sizes, 
(ii) self-gravity in inhomogeneous clumps/cores leading to sub-fragmentation of 
collapsing condensations, (iii) dynamical instabilities in massive disks - formed 
by the initial collapse of rotating cores - fragmenting them into binary or multiple 
star/disk systems, (iv) outward momentum flux from stellar radiation and/or 
MHD winds truncating accretion onto forming stars from the outer parts of 
their parent cores. The relative importance to the final IMF of each of these (and 
other) processes remains to be determined, and major technical challenges are 
involved in attacking any of these questions via direction numerical simulations. 
The large range of scales (nine orders of magnitude!) involved in going from a 
cloud to a star points to the need for adaptive mesh refinement (AMR) schemes 
in order to follow even a portion of the overall process. 

For the present, we can begin by assessing the clumpy structure produced by 
turbulence and self-gravity at moderate scales in uniform-grid MHD simulations 
of GMCs. This clumpy structure can be characterized in many ways, and ap- 
plying varying methods is valuable for understanding and illustrating different 
aspects of the formal dynamical problem, as well as for interpreting observations. 

One class of non-hierarchical structure-analysis methods is similar to the 
CLUMPFIND algorithm introduced by Q. In this method (see ^ for details), 
one chooses a threshold volume density pth (or column density -S'th), identifies 
the set of local maxima in the data cube (potentially first smoothing the data to 
reduce pixel noise) , and then defines clumps by assigning matter at p > pth to the 
nearest local maximum. The shape of a clump can be quantified by computing 
the ratios of principal axes in its moment of inertia tensor, and the importance 
of self-gravity in binding a condensation can be described, e.g., by computing 
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Fig. 3. Number of clumps as a function of mass for two different clump identification 
algorithms. Lower panel shows results of excluding clumps that lie within larger clumps; 
upper panel shows results of counting clumps at any given spatial scale regardless of 
whether they lie within a larger concentration. 



the ratio of gravitational energy to the (weighted) sum of kinetic, thermal, and 
magnetic energies. 

A full discussion of the results of applying this clump-finding method to a set 
of decaying-turbulence (with initial A4 = 14), self-gravitating MHD simulations 
is given in Results from this analysis (taking pth/p = 10) include: 

• At any time, only the high-mass wing of the clump distribution is self- 
gravitating; 

• Characterizing the spectral shape of the high-mass wing as dN/dM oc M~^, 
the slope x is in the range 2-3, becoming shallower over time from mergers; 

• The turnover in the mass spectrum is spatially well-resolved (at ~ 8 grid 
zones across) with Mpcak/Mtot of a few times 10~^; 

• The minimum clump mass is typically a factor 10 below the peak value; 

• Clump shapes are intrinsically triaxial, and project to having two-dimensional 
axis ratios ^2:1; 

• The shapes of the mass functions of apparent clumps in column density maps 
are similar to those of true three-dimensional clumps, but shifted to larger 
masses by an order of magnitude. 

An interesting point is that Afpoak, and also the minimum mass Mmin of a clump 
for given pth/p — 10, do not vary with the value of /3. Also, although the Mach 
number declines by a factor 4-5 over the course of the simulations, the peak 
and minimum masses do not change significantly. This suggests that the mass 
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function of clumps retains a "memory" of the dynamical history of a cloud, 
rather than being determined solely by the cloud's instantaneous turbulent Mach 
number and spectrum (as proposed in 

The results for typical slopes of the high-mass end of the clump mass function 
are intriguingly similar to the value 2.35 for the Salpeter stellar IMF, which 
also appears to describe the core IMF in forming clusters Similar results 
have also been obtained by analysis of a variety of simulations by other groups 
(e.g. [Q, ^^,|4|). While the conclusions from these preliminary analyses are 
promising, many questions still remain open. It not yet clear in general how the 
mass scales (the peak, minimum, and maximum) depend on the input parameters 
{M , as well as the spectral shape, and potentially P) and on dynamical history. 
Uncertainties also remain in how three-dimensional clumps relate to clumps 
seen in projection, and whether the latter distribution can be used to deduce 
the former. 

Another important set of questions is how the definition of a "clump" ~ both 
the specific clump identification algorithm with its chosen set of parameters (such 
as pth), and the overall category of algorithms in which a specific method lies 
- affects the results. A clear categorical distinction is between non-hierarchical 
identification algorithms like CLUMPFIND (in which every mass element is 
assigned to a single clump), and hierarchical algorithms, in which a given mass 
element may be counted as part of many clumps, at different levels of a hierarchy. 
As an astronomical analogy, it is clear that for many purposes it is valuable to 
count galaxies whether or not they are part of larger clusters or supercluster; 
counting "objects" is a function of the spatial scale under consideration. 

As an example of how hierarchical considerations affect the mass spectrum, 
consider the distributions shown in Fig. 3. For this analysis, a "clump" is any 
cubic region at a given spatial scale in which the density exceeds the mass- 
weighted mean density for the ensemble of cubes at that scale (see |^ for 
further details on this "region of contrast" algorithm). For the lower panel in 
Fig. 3, only clumps that do not lie within other clumps are counted for the 
mass spectrum; for the upper panel, the spectrum counts clumps regardless of 
"overlap." As should be expected, allowing for clumps-within-clumps leads to 
a relatively steeper mass spectrum, with a slope dN/dM oc M~^'^ in this in- 
stance. The no-overlap spectrum is dN/dM oc M~^-^. This difference between 
the hierarchical and non-hierarchical mass spectra is very interesting because 
it is reminiscent of the difference between the steep stellar IMF, and shallower 
observationally-determined GMC mass functions (e.g. Q). Because the density 
condensations produced by turbulence as "initial conditions" for collapse are hi- 
erarchically nested, this difference in slopes offers intriguing support for the idea 
that fragmentation during gravitationally-coUapsing stages may play a crucial 
role in defining the stellar IMF. An important question for future AMR simu- 
lations to address is when and why gravity chooses an "inner" versus "outer" 
mass scale for final collapsed objects. 
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4 Linewidth-size Relationships 

The direct observables produced by spectral-line mapping of a molecular cloud 
are data cubes of intensity as a function of two plane-of-sky positions and the 
line-of-sight velocity (in radio astronomy, the intensity is described as a bright- 
ness temperature) . In principal, one would like to extract the spatial distribution 
of velocity and emissivity from these data cubes. Because only projected data 
are available, and a turbulent cloud has no spatial symmetries to exploit, di- 
rect inversion is not possible. However, one may still hope to deduce statistical 
properties of the turbulence from the full intensity data cube. Various complex 
techniques to do this are under development by several groups - including Prin- 
cipal Component Analysis , the Spectral Correlation Function |^ , |^ ; 
and Velocity Channel Analysis P4 , |p3| . Here, I will briefly discuss a simple 
technique to estimate a theoretically fundamental - and observationally much- 
investigated - property of turbulence, the variation of linewidth with physical 
size scale. 

Averaged over volumes with the same size in all three directions, the mean 
linewidth simply reflects the underlying one-dimensional velocity power spec- 
trxim, since 

If the turbulence has a power-law spectrum, w^(fc) cx then iTt,(s) oc 

5(a-3)/2^ so for e.g. Kolmogorov or Burgers spectra with a = 11/3 or 4, a.^ oc 

gl/3 pj. 

In observations, however, any region of size in projection (on the plane of 
the sky) in general extends over a scale at least as large along the line of sight. If 
s is small compared to the overall scale of a cloud, L, then the observed linewidth 
from a region of projected area can have contributions from L/s S> 1 volume 
elements of size along the line of sight. If power increases with scale (a > 3), 
then the velocity centroids of the multiple volume elements on average differ, 
such that the linewidth over the area integrated over the whole line of sight 
will on average exceed (Ti,(s), reaching up to ay{L). However, because of the sta- 
tistical nature of the distribution (e.g. if it obeys Gaussian random statistics), for 
some projected positions the centroids of the multiple s^ volume elements will 
differ very little, such that the line-of-sight integrated velocity dispersion will be 
close to ay{s). From this argument, one expects that the mean linewidth would 
vary weakly with projected size, whereas the lower envelope of the linewidth 
distribution would vary more strongly with projected size, and in fact trace the 
underlying three-dimensional mean linewidth-size relation. Analysis of simula- 
tion data cubes indeed bears out this expectation plf , as shown for example in 
Fig. 4. 

The foregoing discussion is helpful for interpreting well-known observational 
aspects of molecular cloud kinematic scalings. In particular, by relation to Fig. 
4, one may understand why the mean linewidth-size relation for apparent clumps 
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Fig. 4. Linewith-size relations for projected and 3D simulated data. Distribution 
(squares) shows total linewidth for apparent clumps in a map vs. linear dimension 
of projected clump (i.e. square root of its area); dashed line shows (linear) fitted mean 
linewidth-size distribution for projected clump map (slope is 0.2). Solid curve shows 
mean linewidth-size relation for 3D volumes from the same data cube. Data set is 
P = 0.1 model with Mach number 5 from with threshold density p/p — 3. 

observed in molecular tracers with moderate Pcrit/p is relatively flat (e.g. , [p| , p8| ) 
compared to the relatively steeper "Larson's Law" (ay cx s-^^^) linewidth-size re- 
lations (cf. ||l9|, that apply to objects that are spatially "isolated" - either 
because Pcnt/p is large (for dense cores within clouds) or because of phase differ- 
ences with their surroundings (for molecular clouds within the atomic ISM) . For 
tracers with critical density near the (mass-weighted) mean density in a cloud, it 
is not unlikely for multiple structures that are separated and in relative motion 
along the line of sight to contribute to the linewidth at a given position on a map. 
For tracers with higher critical densities, the probability of chance projections of 
dense condensations is much lower because there are many fewer such condensa- 
tions. Spectral correlation function analyses |35j also demonstrate that simulated 
spectral data cubes from higher-density tracers yield more spatial variation than 
do data cubes from lower-density tracers, for the same reason: the spectrum at 
a given point in a map from a low-density tracer samples more completely along 
a line of sight, and is thus more representative of the mean spectrum averaged 
over a whole cloud - compared to the spectrum from a high-density tracer. 

This discussion of linewidth-size relations in 2D and 3D also serves to il- 
lustrate the point that coherent structures in position-velocity space in general 
differ from coherent structures in 3D position (physical) space; detailed analy- 
ses of simulation data cubes have shown this in a number of different ways (e.g. 
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As a consequence, mass functions of density condensations, or other 
measures of structure in the density such as its Fourier power spectrum, can- 
not necessarily be obtained by treating the hne-of-sight velocity as a surrogate 
for line-of-sight position. Instead, it is necessary to use statistical approaches 
to diagnose structure in the physical density distribution - just as statistical 
approaches are required for diagnosing structure in the velocity distribution. 
Depending on the relative power in velocity and density fluctuations, true den- 
sity structure may be best discerned by integrating intensity over velocity and 
then correcting for physical superposition of dense structures using information 
from two-point correlation functions in the column density (see [|22| for related 
discussion). 

5 Polarization as a Magnetic Field Diagnostic 

It has long been recognized that polarization studies are important for diagnos- 
ing basic properties and structure of the ISM, because they provide relatively 
direct access to the elusive - but dynamically important - magnetic field (e.g. 
[ p^ , 1^ ) . Provided that dust grains preferentially align with short axes parallel 
to the local direction of the magnetic field (|l^; see e.g. recent review of |^ for 
thorough discussion of alignment mechanisms), an ordered B-field will lead to 
observable polarization. For grains aligned with their short axes parallel to the lo- 
cal magnetic field, background stars are observed in optical/near-IR wavelengths 
with polarization parallel to the local magnetic field, and local dust emission is 
observed in far-IR/sub-mm wavelengths with polarization perpendicular to the 
local magnetic field. If the angle of B with respect to the plane of the sky is 
i, the local contribution to polarization is cx cos'^ i times the difference between 
long- and short-axis grain crossections times the density of polarizing grains. 
Taking the simplest-possible assumption of a uniform ratio of polarizing-grain 
density to gas density (but see below), it is straightforward to create simulated 
polarization maps by integrating the radiative transfer equations for the Stokes 
parameters ( |^ , ||2^ , j24| ) . For a medium with optical depth r ^ 1, the fractional 
polarization in (thermal) emission and in dust-absorbed starlight are related by 
-Pom — Pahs/T (e.g. ||l7|| ), SO that polarization from emission is proportional to 
polarization from extinction divided by the column of intervening matter. 

Maps and analyses of polarized extinction and emission [^,||lj| com- 
puted from MHD simulations have recently been presented by several groups. 
One question of interest is how the fractional polarization varies with the col- 
umn N = J nds of absorbing or emitting material. For a uniform magnetic field, 
uniform polarization efficiency of grains, but spatially-nonuniform distribution 
of matter, Pabs would increase linearly with N , while Pom would be independent 
of N . For a spatially nonuniform field, the situation is much more complicated. 
In this case, the increase in Pabs or Pcm/T with N would (a) in general be shal- 
lower because variations in the magnetic field direction decorrelate the direction 
of grains, reducing the net contribution to polarization per unit length along 
the line of sight, and (b) no longer follow a linear relation if the amplitude of 
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Fig. 5. Distributions of simulated polarized-extinction as a function of column density, 
for decaying-turbulence simulations with A4=7 from [^. Scatter plots show Pabs for a 
randomly-chosen subset of positions on the projected map, for three different projection 
directions, and three different mean magnetic field strengths. Squares with 1-cr error 
bars show the mean P — N relation for the full data sets. Dotted lines show what the 
P — N relation would be for uniform magnetic fields. P is arbitrarily normalized to a 
value 0.1 for uniform density N — N for a projection perpendicular to Bo 



field fluctuations is large and/or if the number of effective correlation lengths 
of magnetic field orientation along the line of sight varies with column density 
(this can yield, e.g., Pabs « N^^'^ fUl)- For weaker mean magnetic fields Bq 
and a given power spectrum of fluctuations 6B, the directional decorrelation in 
polarization occurs at a physically smaller scale than for the strong- i?o case, so 
that lower polarization is expected for a given column density. These expected 
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trends are indeed evident in distributions of Pabs vs. column density from MHD 
simulations, as shown for example in Fig. 5. 

Interestingly, observed distributions (e.g.[|l2|, j^) of polarized-extinction vs. 
column (or Ay) in molecular clouds do not show the behavior evident in Fig. 
5, which would correspond to a secular increase in Pabs up to Ay of 30 or more 



(for N corresponding to a typical GMC Ay ~ 7.5 l27|). Instead, the observed 
increase of Pabs with Ay flattens near Ay of unity, possibly indicating that grain 
alignment fails in the deep interiors of clouds (e.g. [pOl). Additional support for 



this idea comes from comparison of observed distributions (e.g. |15|) of Pem 
vs. intensity (proportional to N) in dense cores with simulated distributions; 
the simulations show insufficient decrease in Pem with column unless high- Ay 



regions have decreased polarizing efficiency [36 



Potentially, one of the most important applications for polarization studies 
is to use the variation in polarization directions to diagnose the strength of the 
magnetic field |^ . The basic physical idea behind the so-called "Chandrasekhar- 
Fermi" method is that weaker magnetic fields produce lower tension forces for 
a given displacement, so that for a given velocity field, a lower mean magnetic 
field strength results in larger distortions in the magnetic field direction. For 
the case of a single low-amplitude Alfven wave in a uniform magnetic field with 
plane-of-sky component Pp, the magnetic field and velocity perturbations obey 
SB /Bp = Sv/vA,p, where ua,p = Pp(47rp)~^/^. For a wave with perturbation 
direction in the plane of the sky, the dispersion in polarization directions is 
S(j) = {{5B / BpY)^/"^ if the net polarization is either parallel or perpendicular to 
B. If there is an Alfven wave component with perturbations along the line of 
sight having the same amplitude as the component with perturbations in the 
plane of the sky, then Sv/vA,p = (^ios)^^^/''^A,p — Scj), so that Bp is given by 
PcF = {inpY^^Sv/ScI). Thus, with measures of the mean density, the observed 
velocity dispersion, and the dispersion in polarization angles, the plane-of-sky 
field strength may in principle be estimated. 

Since several of the idealizations described above do not hold in a real molec- 
ular cloud, it is desirable to test and/or recalibrate the Chandrasekhar- Fermi 
(hereafter "C-F" ) relation using MHD simulations, which allow for more complex 
dynamical structure. Figure 6 shows an example of such a test using simulated 
polarized-extinction data, comparing the true mean plane-of-sky field strength 
with the "one-wave/equipartition" C-F estimate given above. As reported in 
when the dispersion in polarization angles is sufficiently small {S(f> < 25°) 

- such that linear theory is adequate, a good estimate of the plane-of-sky field 
is (Bp) ~ 0.5PcF- Evidently, the presence of more than one wave along the line 
of sight reduces (50, so that Bcf tends to overestimate the true Bp. For larger 
dispersions in polarization angle. Fig. 6 shows that the measure O.bBcp "cali- 
brated" at small Scf) may either under- or over-estimate Bp. Note that even if the 
mean magnetic field is strong, Scf) can be large - and linear theory inappropriate 

- for system orientations in which the mean magnetic field direction lies near 
the line of sight. Analyses of simulated polarized-emission maps, combined with 
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Fig. 6. Plane-of-slcy component of the mean magnetic field, 
Chandrasekhar- Fermi estimate Bcf, for A4=7 (four-sided symbols) and A4—9 (five- 
sided symbols) decaying turbulence MHD simulations |^]; 5(j) is the corresponding 
dispersion (in degrees) of directions in the polarized-extinction map. Solid, starred, 
and open symbols represent j3 = 0.01, 0.1, 1.0 models. 



kinematic measurements |lj] , or with synthetic radiative transfer spectral maps 
, yield similar results for testing of the C-F formula. 

We note that the C-F method measures only the mean plane- of-sky mag- 
netic field (Bp)] Figure 7(a) shows that even for small dispersion in polarization 
angles, there can be large variation in the total field strength compared to the 
"corrected" C-F estimate. In any individual cloud, the line-of-sight field can be 
estimated by the Zeeman effect (at least in principle; in practice this is diffi- 
cult), with the two estimates combined to obtain the total three-dimensional 
field strength (e.g. Another potential caveat for observational application 

of the C-F method is that the principal correlation scale for the dispersion in 
magnetic field directions must be resolved in the plane of sky by the effective 
observational "beam" . For polarized extinction, this poses no difficulties be- 
cause background star observations have minimal beam thickness. For polarized 
emission from warm cores, however, large-beam averaging of spatial fluctuations 
tends to reduce 6(j) relative to its value for a "pencil-beam" observation, poten- 
tially resulting in an overestimate of Bp unless an appropriate correction factor 
is apphed [l|]. 

When variations in the polarization directions in a map are large, the fluctua- 
tions of the plane-of-sky magnetic field must be comparable to the mean value of 
Bp. In this case, it is not possible to distinguish from polarization studies alone 
whether (a) a cloud has a large mean magnetic field that is "hidden" along the 



Cloud Diagnostics from MHD Simulations 15 



]— I— I— I— P 



1.5 



0.5 



^ 



4 

O 



I 



1 ^ 


- (b) - 

■ 


1.5 


1 

# 






CO 


i 


1 




_ >* 


□ 


0.5 


□ 








20 40 

60(deg) 



60 



Fig. 7. (a) Total (rms) magnetic field strength relative to C-F estimate Bcf vs. disper- 
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line of sight (as is the case for the solid points in Fig. 6 having large S(j)), or (b) 
the mean magnetic field is simply weak (as is the case for the starred and open 
points in Fig. 6). Observed line-of-sight velocities combined with an assumption 
of equipartition between magnetic and kinetic energies can still yield an approxi- 
mate measure of the rms magnetic field strength (see e.g. Fig. 7b) that is correct 
within a factor 2, although an important caveat is that this estimate could be 
arbitrarily far off for cases with very strong mean magnetic fields (ua ^ Vrms) 
that happen to lie near the line of sight. [] Fitting formulae that interpolate 
between the small-(50 and laTge-dcj) limits (with the implicit assumption that wa 
along the line of sight is not very large) have been provided by . 



Also note that while Brms is a dynamically-important quantity, it is not equivalent 
to the mean magnetic field strength that enters into the mass-to-flux ratio, 

which ultimately determines whether a cloud or core is super- or sub-critical. 
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6 Summary 

To understand the intrinsic nature of MHD turbulence and the relation between 
the turbulence observed in astronomical systems and the turbulence simulated 
numerically, it is crucial to develop structural diagnostics. These diagnostics 
may: 

• Enhance conceptual understanding of the turbulent phenomenon; 

• Enable determination of astronomical systems' properties that are either 
difficult to observe at all (e.g. B) or only indirectly observable because of 
projection effects (e.g. Vk, Pk)] 

• Provide a physical basis or interpretation for observed empirical "laws" 
(e.g. column density distributions, linewidth-size relations, mass functions 
of clouds, clumps, and stars); 

• Indentify when additional physical ingredients may be needed in a compu- 
tational model. 

From the examples outlined in this paper, it is clear that significant advances 
along these lines have already been accomplished. 

Conceptual and (approximate) quantitative understanding of what deter- 
mines the PDFs of density and column density observed in molecular clouds 
have already been obtained from analyses to date (see §2). The development 
of lognormal statistics from Gaussian random processes under near-isothermal 
conditions appears to obtain robustly for a variety of conditions, with little sen- 
sitivity to the magnetic field strength. Because one-point density statistics are 
subject to "cosmic variance" (if low wavenumbers dominate the power spectrum, 
different realizations of a given power spectrum have significant variation), and 
because degeneracies prevent the direct inversion of column density statistics 
to obtain volume density statistics, there may not be a highly accurate way 
to determine the physically-important volume-averaged density in a cloud from 
more direct observables such as the observed velocity dispersion and observed 
distribution of extinction. Future work that combines one-point column density 
statistics with two-point correlation functions could potentially be valuable in 
constraining the spatial power spectrum of density fluctuations. Such analyses 
could also be useful in relating mass functions of clumps seen in projection with 
true three-dimensional clumps (see §3, §4). 

Analysis of clumpy structure in model clouds (see §3) shows intriguing corre- 
spondence to observations: mass functions of self-gravitating structures are com- 
parable to the Salpeter IMF, and there are significant differences between steeper 
(more "stellar- like" ) mass functions when subcondensations are not subsumed 
into larger structures, and shallower (more "cloud-like") mass functions when 
substructure is discounted. These findings support the concepts that clumping 
imposed by turbulence, as well as fragmentation during gravitational collapse, 
may both be important in establishing the IMF. Future work is required to de- 
termine whether there might be a relationship between Mach numbers (instanta- 
neous and historical), the overall size and mass of a cloud, and the characteristic 
"peak" sizes and masses of self-gravitating condensations within a simulation. 
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More coverage of model parameter space (allowing for different power spectra, 
forcing, etc.) will be important for deciding how sensitive the resulting clump 
mass functions may be to physical conditions. 

Some of the longest-established empirical results about molecular clouds con- 
cern the correlations among physical scale and spectral linewidth, and analyses 
of simulations have proven valuable in interpreting how these empirical "laws" 
relate to the underlying properties of turbulent clouds and how they are observed 
(sec §4). Steep ("Larson's law") dependence of linewidth on size probably reflects 
the true three-dimensional power spectrum, since the structures to which these 
steep laws apply are observed in tracers that render them well-separated from 
the background. On the other hand, the weak dependence of mean linewidth 
on 2D size of apparent moderate-density clumps within clouds may largely arise 
from projection effects - with the scale sampled by the velocity dispersion on 
average extending over much of the whole cloud's depth. Both simple methods 
using the lower envelope of the linewidth-size distribution - and more compli- 
cated methods using detailed spectral shapes and their spatial correlations, are 
very promising for being able to distinguish the true three-dimensional power 
spectrum from molcciilar line data cubes. 

Since magnetic fields are difficult to measure directly, there are particularly 
strong motivations to develop indirect diagnostics (see §5). Polarization either 
in absorption or emission is sensitive to the local direction of the magnetic field, 
and the variation in the local direction of the magnetic field is sensitive to the 
magnitude of the magnetic field and level of turbulence in a cloud. Thus, one 
might hope to combine observed measures of variation in the polarization direc- 
tion with observed measures of turbulence via molecular linewidths to infer the 
mean magnetic field strength; this kind of indirect method was originally pro- 
posed by Chandrasckhar and Fermi. Testing the Chandrasckhar-Fcrmi method 
with simulation "data" shows that for low dispersion in the polarization angle, 
recalibration by a factor one- half indeed yields a good measure of the mean plane- 
of-sky magnetic field. When the polarization direction has large fluctuations, the 
Chandrasekhar-Fermi method loses accuracy, but simulations show that an as- 
sumption of magnetic/kinetic equipartition is usually correct within a factor 
^ 2. Because undcrresolution tends to enhance the estimated field strength, and 
because high-column/high-density regions may not be efficient polarizers, there 
are some potential caveats in applying the Chandrasekhar-Fermi method for 
polarized-emission data. The promise shown by analyses to date, together with 
the potential to obtain large-scale maps of polarized absorption and emission in 
molecular clouds, marks this area as an important direction for further research. 

These detailed results are exciting, and represent only a small sample of the 
progress that has taken place in this field to date. Perhaps the most fundamental 
advance, however, has been the movement to match our sophisticated concept 
of what a molecular cloud is - a, complex structure with multiple-scale, large- 
amplitude, turbulent fluctuations in all fluid variables; with a commensurately 
sophisticated way to diagnose structure - by developing and testing analytical 
tools on detailed simulation data cubes that are self-consistent, time-dependent 
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solutions of the MHD equations. Recent accomplishments have greatly advanced 
our field, but we are still very much in the era of discovery. Ongoing interplay 
between simulation, analysis, and observation will be essential for continued 
progress in building a comprehensive dynamical model of molecular clouds. 



I am grateful to J. Stone and C. Gammic for permission to present results 
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